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ABSTRACT 

We investigate the origin of the high Faraday rotation measures (RMs) found for polarized radio 
galaxies in clusters. The two most likely origins are, magnetic fields local to the source, or magnetic 
fields located in the foreground intra-cluster medium (ICM). The latter is identified as the null hypothesis. 
Rudnick & Blundell (2003) have recently suggested that the presence of magnetic fields local to the source 
may be revealed in correlations of the position angles (PAs) of the source intrinsic linear polarization 
and the RMs. We investigate the claim of Rudnick & Blundell to have found a relationship between the 
intrinsic PAg of the radio source PKS 1246-410 and its RM, by testing the clustering strength of the 
PAq-RM scatter plot. We show that the claimed relationship is an artifact of an improperly performed 
null-experiment. We describe a gradient alignment statistic aimed at finding correlations between PAq 
and RM maps. This statistic does not require any null-experiment since it gives a unique (zero) result 
in the case of uncorrelated maps. We apply it to a number of extended radio sources in galaxy clusters 
(PKS 1246-410, Cygnus A, Hydra A, and 3C465). In no case is a significant large-scale alignment of PAq 
and RM maps detected. We find significant small-scale co-alignment in all cases, but we are able to fully 
identify this with map making artifacts through a suitable statistical test. We conclude that there is 
presently no existing evidence for Faraday rotation local to radio lobes. Given the existing independent 
pieces of evidence, we favor the null hypothesis that the observed Faraday screens are produced by 
intracluster magnetic fields. 

Subject headings: galaxies: clusters: general — magnetic fields — polarization — radio continuum: 
galaxies 



1. INTRODUCTION 

A detailed summary of the observational evidence for 
the presence of magnetic fields embedded within the ther- 
mal gas in clusters of galaxies is presented in two recent 
review papers (CariUi & Taylor 2002, and Widrow 2002) 
and is only briefly summarized here. 

Magnetic fields in galaxy clusters are known to exist due 
to the detection of diffuse cluster wide synchrotron emis- 
sion (Willson 1970) in a number of clusters. This emission 
is not associated with individual galaxies in the cluster 
and is observationally classified as halo or relic emission 
(Feretti 1999). Although initially thought to be relatively 
rare objects, recent surveys have significantly increased the 
number of radio halos and relics (Giovannini, Tordi, & Fer- 
etti 1999; Giovannini & Feretti 2000; Kempner & Sarazin 
2001; Bacchi et al. 2003). The presence of this diffuse 
synchrotron emission reveals the widespread distribution 
of magnetic fields within the intracluster medium (ICM) 
in these clusters. Equipartition assumptions provide min- 
imum energy magnetic field estimates in radio halos of 0.1 
- 1 /iC (Feretti 1999; Bacchi et al. 2003) and 0.4 - 2.7 fxG 
in radio relics (Enfilin et al. 1998). 

Estimates of volume averaged intracluster magnetic field 
strengths can be obtained by comparing the synchrotron 
and inverse Compton emission (e.g., Harris & Grindlay 
1979, Rephaeli et al. 1987). The same relativistic particle 
population which produces the diffuse synchrotron emis- 
sion is expected to up-scatter the ambient photon field in 
the ICM to produce inverse Compton X-rays. Typical esti- 



mates yield a volume-averaged intracluster magnetic field 
in the range of 0.2 - 1 /iC (see e.g. CarilH & Taylor 2002 
and references therein) . We emphasize that inverse Comp- 
ton based estimates have to be regarded as lower limits. 

Faraday rotation measure (RM) studies of extended ra- 
dio sources embedded in clusters probe the magnetic field 
component integrated along the line-of-sight, weighted by 
the plasma electron density. Faraday rotation is the ro- 
tation of the position angle (PA) of the plane of linear 
polarization of radio waves, which traverse a magnetized 
plasma. The RM is the proportionality constant of this 
A^-dependent effect: PA(A) = PAq + RM, where PAq is 
the source intrinsic PA, only directly observable at short- 
est wavelength. The RM and PAo can be obtained by 
multi-frequency measurements of a polarized radio source. 

CariUi & Taylor (2002) review the generally large RM 
values obtained from extended radio sources embedded in 
clusters, and summarize the evidence that indicates that 
the observed RM of embedded sources are indeed probing 
a foreground ICM, with conservative estimates of mag- 
netic field strengths between a few fiG and 10s of ^G. 
Similar magnetic field strengths are determined from sta- 
tistical RM studies by Kim, Tribble k Kronberg (1991), 
and Clarke, Kronberg & Bohringer (2001). 

Finally, the asymmetric depolarization of double radio 
lobes embedded in galaxy clusters can be understood as 
resulting from a difference in the Faraday depth of the two 
lobes (Laing 1988; Garrington et al. 1988) which strongly 
supports the association of RMs with ICM magnetic fields. 
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We note that such a difference in Faraday depth might 
also be explained by local effects near the source, if there 
is a difference in the mixing layers of the magnetized radio 
plasma with the dense environmental gas between the ra- 
dio lobe head side and the back-flow side of an FR 11 radio 
galaxy. In contrast, this scenario would have difficulties to 
explain the observed strong RM and depolarization asym- 
metry of FR I radio galaxies like Hydra A (Taylor & Perley 
1993), which are not believed to have back-flows. However, 
this discussion shows that the association of the RM with 
the intra-cluster medium - in the following regarded as 
the null hypothesis - is not unambiguous, since it could 
also be produced in a magnetized plasma skin or mixing 
layer of the observed radio galaxy (Bicknell et al. 1990). 
For that reason Rudnick & Blundell (2003, hereafter RB) 
used Faraday rotation measure observations of PKS 1246- 
410 by Taylor, Fabian & Allen (2002) in an attempt to 
estimate the fraction of the measured RM signal which is 
local to the radio source. If most of the RM signal is local 
to the radio source then the derived ICM magnetic field 
strength should be significantly lower than if all of the RM 
originates in the ICM. 

In the following we examine RB's claim for evidence 
of source-local magnetic fields. We do this first through 
a general argument (Sect. 2) followed by an illustrative 
numerical example (Sect. 3). Then, in order to have an 
unbiased tool to measure cross correlations of RM and PA 
maps, we specify the mathematical requirements of a suit- 
able statistic (Sect. 4), and construct a gradient alignment 
statistic which meets our quality criteria (Sect. 5). This 
statistic is also sensitive to correlated noise resulting from 
the map making imperfections (Appendix B). We there- 
fore design a gradient vector product statistic, which is only 
sensitive to such imperfections, in order to separate spuri- 
ous signals from any astrophysical signal (Sect. 6). Appli- 
cation of our approach to several datasets (PKS 1246-410, 
Cygnus A, Hydra A, and 3C465) does not reveal any evi- 
dence for source-local magnetic fields (Sect. 7). We con- 
clude (Sect. 8) that the null hypothesis that the RMs are 
connected to the ICM is not only consistent with, but also 
strongly favored by present data. The basic assumption 
for Faraday rotation based ICM magnetic field estimates 
of RM being mostly generated by the ICM magnetic fields 
seems therefore to be valid. 

2. PAo-RM SCATTER PLOTS 

If the Faraday rotation observed in cluster radio galaxies 
is produced by a mixing layer between the radio lobe and 
the ICM gas then it is possible that the magnetic struc- 
tures within the Faraday screen are somehow related to 
the magnetic field orientation within the radio lobes. In 
such a case, co-spatial structures in RM and PAq maps are 
possible since both quantities contain geometrical informa- 
tion about the source local magnetic field geometry: The 
PAo gives the plane-of-the-sky direction perpendicular to 
the sky-projected magnetic field within the radio lobe, and 
the RM the line-of-sight component of the magnetic field 
within the Faraday rotating medium. 

The idea of RB is to test the above scenario by searching 
for the expected co-spatial structures in the PAg and RM 
maps of the radio galaxy PKS 1246-410. One difficulty 
with such an analysis is that no direct correlation between 



the PAq and RM values can be expected even for co-spatial 
structures, since there is no generic reference point for 
PAqs. Therefore RB analyzed the PAq-RM scatter plot 
generated by plotting for each map pixel the PAq and the 
RM value in the same diagram. They argue that local co- 
alignment of PAq and RM values should lead to a strongly 
clustered point distribution in the scatter plot. Since sta- 
tistically independent PAq and RM maps may also pro- 
duce such clusters they perform a null-experiment. They 
generate synthetic RM maps, which have the same power- 
spectrum as the observed RM map, but random phases. 
With these simulated maps they repeat their analysis and 
find much less clustering in the scatter plot. From this 
they conclude that there is a strong co-alignment in the 
PAo-RM maps of PKS 1246-410, and therefore most of the 
RM should be local to the radio source. In order to give 
further support to their claim, they report a number of re- 
gions in a larger sample of radio galaxy maps where they 
see corresponding structures in the RM and PAq maps. 

This latter finding is statistically very questionable, 
since within any two large datasets the human eye of- 
ten finds apparent correlations within some sub-regions, 
even if the two datasets are uncorrelated. Furthermore, 
observed RM and PAq always carry correlated noise, as 
recognized by RB and analytically demonstrated in Ap- 
pendix B, since they are generated from the same set of PA 
maps. In the worst case the noise produces step-function 
like artifacts at the same location in RM and PAo maps. 
We note that correlated steps are present within the afore- 
mentioned regions of RM and PAq maps and give - at least 
to our eyes the dominant contribution to the PAq-RM 
correlation impression reported by RB. 

However, the claim of a statistical detection of a PAo- 
RM correlation can be rigorously investigated. The null- 
experiment applied by RB is designed to have the same 
two-point autocorrelation function as the observational 
data, but higher order correlations are neglected by the 
random phase realization. Unfortunately, the chosen 
statistic is very sensitive to such higher order correlations. 
Any clustering in the PAq-RM scatter plot is caused by 
having patches of nearly constant values in the PAq and 
RM maps. The appearance of patches requires special re- 
lations of the Fourier phases, or - equivalently - non-trivial 
higher order correlation functions to be specified. 

The strength of the clustering in the scatter plot is in- 
deed strongest if the PAq and RM structures are corre- 
lated. However, the clustering does not disappear if the 
RM and PAq patches have independent distributions, since 
every PAq patch is still overlaid by a small number of RM 
patches, so that the associated clustering in the scatter 
plot only gets split into a corresponding number of smaller 
clusters. These clusters happen to be co-aligned on a ver- 
tical constant PAq line in the plot, since they all belong 
to the same PAq coherence patch. A corresponding effect 
splits the pixels of an RM cell into a horizontal line of 
clusters. Such horizontal and vertical structures are in- 
deed visible in the scatter plot of PKS 1246-410 (Fig. Bl, 
also Fig. 4 of RB). 

We note that a proper null-experiment, which would 
maintain the higher-order correlations, can also be con- 
structed empirically by exchanging subregions of the ob- 
served RM image from one radio lobe to the other. This 
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should keep the same RM correlation functions, but will 
destroy any real PAq-RM alignment since the different 
regions of the source are independent. We have per- 
formed this lobe switching by simply dividing the source 
in roughly two equal regions about the center and shift- 
ing the coordinates in right ascension such that the two 
subregions overlap. We have then plotted the PAq of the 
east (west) lobe against the RM of the west (cast) lobe. 
This source division is preferable to a reflection about the 
declination axis through the core as our lobe shifting re- 
sults in a correlation of central source region PAq with the 
outer lobe region of RM thus eliminating possible radial 
influences on the correlation. The results of such an exper- 
iment are also presented in Fig. Bl. Clustering of points 
in the PAq-RM scatter plots are clearly seen in both plots 
in contrast to RB's claim that clustering results from a 
relation between PAq and RM. 

While RB noted that higher order correlations in inde- 
pendent RM and PAq maps can systematically produce 
spurious signals in their statistic. It seems to follow that 
this implies that their statistic does not allow any conclu- 
sions about the existence of PAq-RM co-alignment. They 
argue that the presentation of a counter example, as we 
provide in the next section, does not invalidate their ap- 
proach, since it cannot be demonstrated that the example 
used is realized in nature. However, we note that RBs 
approach will not be useful if patchiness is important. Al- 
ready an inspection by eye of the maps of PKS 1246-410 
in comparison to RB's random phase maps shows that na- 
ture (or the map making software) produces patchy maps 
which must exhibit clustering in PAq-RM scatter plots, 
independent of whether PAq-RM alignment is present or 
not. 

In the next section, we provide a patchy map generating 
algorithm, which illustrates our general argument. More 

importantly, it allows us to test the gradient alignment 
statistic, introduced in Sect. 5. 

3. SYNTHETIC PATCHY MAPS 

For illustration purposes we construct a simple exam- 
ple, which captures the main effects, but is not meant to 
be exact in all respects. We construct patchy RM and PAq 
maps, which are statistically independent from each other. 
For both maps we use a similar recipe*. 

First a number N of random seed points Xi {i G 
{l...A^}) within a square area is drawn, and then the area 
is split into cells around the seed points by means of a 
Voronoi-tessellation: Each point x of the area belongs to 
the cell of its nearest seed Xi. Then each seed is attributed 
a random value ipi {tjj stands in the following for both RM 
and PAq) and a small two-dimensional random vector Vi/^i 
(= auxiliary RM or PAq gradients within the patches, only 
used for the map construction). Each pixel x within the 
cell of seed i gets a value 



tp{x) = Vi + VtZ-i • (f - Xi) + a{x) , 



(1) 



where a{x) is a small random noise term. The resulting 
map consists of patches with nearly constant values, but 
which exhibit some internal trends and noise. The adopted 
parameters are described in Appendix A. 

^ The computer code can be requested from TAE (ensslinOmpa-garching, 



The PAq and RM maps are slightly smoothed, a 20% 
border region is cut away in order to suppress edge effects, 
and a PAq-RM scatter plot is generated. A typical realiza- 
tion of such a scatter plot is shown in Fig. Bl. A strongly 
clustered distribution is visible even though the individual 
PAq and RM maps were completely independent. Further- 
more, nearly horizontal and vertical chains of clusters are 
visible for the reason given in Sect. 2. The map smooth- 
ing produces bridges between these clusters, since it gives 
intermediate values to pixels which are at boundaries of 
PAq or RM ceUs. 

The deviations from the strict horizontal and vertical di- 
rections visible in the scatter plot of PKS 1246-410 should 
be caused by trends within the coherence cells. We note 
that such structures arc also visible in the simulated scat- 
ter plots of RB, although the smooth realizations of their 
RM maps has smeared them out (see their Fig. 2 for a 
comparison of the patchiness of observed and simulated 
RM maps). 

We therefore conclude that the data of PKS 1246-410 
favors statistically independent PAq and RM maps. The 
occurrence of vertical and horizontal lines in the scatter 
plot of the observational data demonstrates that the PAq 
and RM patches are indeed misaligned. 

In order to further investigate this, we construct a model 
in which the PAq and RM patch positions are absolutely 
identical. We construct this by moving the Voronoi- 
tessellation secid points of the RM map to the locations 
of nearby seed points of the PAq map, thereby assuring 
that there is a one-to-one mapping. All other variables 
(RMi, PAo i etc) were kept as before. The recomputed 
RM map has therefore exactly the same patch locations 
as the PAq map. The horizontal and vertical cluster align- 
ments and stripes are absent there (see Fig. Bl). There 
are now several stripes with diagonal orientations due to 
pixels at the PAq-RM cell boundaries which received si- 
multaneously intermediate values in PAq and RM by the 
smoothing. 

4. REQUIREMENTS FOR SUITABLE STATISTICS 

In order to quantitatively investigate the location of the 
Faraday rotation measure screen we must develop a statis- 
tical basis for distinguishing local and foreground effects. 
Therefore, we have to construct a statistic which allows 
us to search for local correlations between PAq and RM 
without requiring a global relation between these quanti- 
ties. In order to design a proper statistic one should first 
specify mathematical requirements in order to avoid po- 
tential pitfalls. We require any suitable statistic A to fulfill 
the following conditions: 

1. The statistic should be sensitive to the presence 
of correlated spatial changes of PAq and RM, in- 
dependent of the local values of these quantities. 
Correlations in PAq and RM maps can not be ex- 
pected to follow a functional form like RM(af) = 
f{PAo{x)) + 6f{x). 

2. The statistic A should not require a null experi- 
ment, which is always problematic since the con- 
struction of a synthetic dataset with all the impor- 
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tant statistical properties 
duced is very difficult. 



of the real data repro- 



3. The statistic has to provide unique expectation val- 
ues in case of uncorrelated {A = 0) and in case of 
fully correlated {A = 1) maps. 

4. The statistic should be analytic and sufhciently sim- 
ple, allowing that basic properties can be derived 
and understood analytically. 

5. The significance of the statistic should increase 
monotonically with the map size. 

The approach of RB to calculate the autocorrelation of 
PAq-RM scatter plots fulfills our requirement No. 1, or 
- more exactly - requirement No. 1 is inspired by their 
work. However, their approach requires a null experiment 
(violating No. 2), since it is not clear a-priori what the 
meaning of the derived correlation in the scatter plot is 
(violating No. 3). Although it is possible to derive analytic 
equations describing the autocorrelation of a scatter plot of 
two quantities, RB's application of a median weight filter 
to the RM and PA maps is a strongly non-linear operation 
which makes it very hard to understand the mathemati- 
cal properties of their data treatment (violating No. 4). 
Finally, if one imagines the map size growing to infinity, 
one easily realizes that the scatter plot will be smoothly 
filled with points, independent if there were local correla- 
tions between PAq-RM or not. Therefore the significance 
of the method of RB starts to decrease when the clumps 
are so densely spaced that they start to merge, and van- 
ishes in the limit of an infinitely large map (violating No. 
5). Applying a median weight filter to maps as RB did 
sharpens the clumps in the scatter plot considerably, thus 
suppresses the clump merging and therefore enhances the 
autocorrelation signal. 

5. GRADIENT ALIGNMENT STATISTICS 

In order to construct a statistic which fulfills our re- 
quirements, we introduce the gradient alignment statistic 
A of different maps by comparing the gradient p = V RM 
and q = VPAq.^ The idea is to check for alignment of p 
and q, indicating correlated changes in PAg and RM, thus 
fulfilling condition No. 1. Since the absolute values of RM 
and PAq are not of any significance for the question of co- 
alignment, the comparison should give the same signal for 
parallel and anti-parallel gradients. Therefore, instead of 
the scalar product 

p-q = pq (cos (j)p cos 4>q -|-sin (j)p sin = p^ q^ +py Qy (2) 

of p = {px;Py) = p(cos(/)p,sin^p) and q = {qx,qy) = 
q (cos (f>q, sin (/)q) we construct an alignment product 

(PjQ) = PQ (cos 2(f)p cos 2(f>q + sin 2(/)p sin 2(j)q) 



The alignment product has the properties (j),ap) = \a\p'^ 
for a being any real (positive or negative) number, and 
(PtO) = —pq a p -i- q- An isotropic average of the align- 
ment product of two 2-dimensional vectors leads to a zero 
signal, since the positive (aligned) and negative (orthogo- 
nal) contributions cancel each other exactly. Therefore our 
requirement No. 2 for suitable statistics can be guaranteed 
by the usage of the alignment product. 

We then define the alignment statistics of the vector 
fields p{x) = VRM(f) and q = VPAo(f) as 



A = A[p,^ 



jd^x\m\\m\ 



(4) 



which has all required properties: 

No. 1: The statistic does not depend on any global 
relation between PAq and RM. This can be demonstrated 
by changing a potential functional dependence using non- 
linear data transformations. Any non-pathological^, piece- 
wise continuous and piecewise monotonic pair of transfor- 
mations RM* = ^(RM) and PA^ = T(PAo) do not destroy 
the alignment signal due to the identity 

(VRM*, VPAo*) = |S"(RM)T'(PAo)| (VRM, VPAq). 

(5) 

Inserted in Eq. 4 one finds that the weights of the different 

contributions to the alignment signal might be changed by 
the transformation, but except for pathological cases any 
existing alignment signal survives the transformation and 
no spurious signal is produced in the case of the null hy- 
pothesis of uncorrelated maps. 

No. 2 & 3: A simple calculation shows that the ex- 
pectation value for A is zero for independent maps and it 
is unity for aligned maps. For illustration, the simulated 
pair of independent maps has A = —0.03, which can be 
regarded as a null-experiment testing our statistic with a 
case where RB's statistic incorrectly detects co-alignment. 
The simulated pair of co-aligned maps has A = 0.89 illus- 
trating ^'s ability to detect correlations. Property No. 3 
implies requirement No. 2. 

No. 4: From the discussion so far it should be obvious 
that many essential properties of the alignment statistics 
can be derived analytically. However, as an additional 
useful example we estimate the effect of a small amount 
of noise present in both maps. Each noise component is 
assumed to be uncorrelated with RM, to PAq, and also to 
the other noise component. We find for small noise levels 



..^ ^ , „^ A[p,q\ 

A\p + Sp,q + 5q\^ . . 



(1 + (5pVp2) (1 + Sq2/q2^ 



(6) 



{pI - pl){ql - ql) + ^PxPy qx qy 



where the bar denotes the statistical average. This re- 
lation holds only approximatively, since the non-linearity 
of A prevents exact estimates without specifying the full 
probability distribution of the fluctuations. As can be seen 
from Eq. 6, uncorrelated noise reduces the alignment sig- 
(3) nal, but does not produce a spurious alignment signal, in 
contrast to correlated noise, which usually does. 

^ We calculate gradients of a quantity i/j defined on our maps by assigning the pixel position with S/ili = {ipi^ij + — V'lJ ~ 

ipij^i + ~ ~ "01+1, j): where Vij denotes the value of ^p at the pixel position Neither the small diagonal shift by 

1/2 pixel in i and j directions, nor the missing normalization of the so defined gradient have any effect on our statistics. PAo-gradients are 
calculated using subtraction modulo 180° in order to account for the cyclic nature of PAs. We note that we apply this scheme to the maps 
directly, even though for our synthetic maps gradient-like auxiliary quantities were defined and used during their construction. 
® A pathological transformation would e.g. split the RM or PA value range into tiny intervals, and randomly exchange them or map them all 
onto the same interval. 
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No. 5: If the RM and PAq maps are enlarged by an 
additional region with the same statistical properties the 
expectation value of A is unchanged, due to the averaging 
property of Eq. 4, and fluctuations in A decrease due to 
the central limit theorem. 

6. DETECTING MAP-MAKING ARTIFACTS 

Before the gradient alignment statistic can be applied 
to real datasets, it has to be noted that it is very sensitive 
to correlations on small scales, since it is a gradient square 
statistic. Observed RM and PAq maps will always have 
some correlated fluctuations on small scales, since they 
are both derived from the same set of radio maps, so that 
any imperfection in the map making process leads to cor- 
related fluctuations in both maps. Therefore the strength 
of signal contamination by such correlated errors has to be 
estimated, and - if possible - minimized before any reliable 
statement about possible intrinsic correlations of PAg and 
RM can be made. This point can not be overemphasized! 

The level of expected correlated noise is investigated in 
Appendix B. Fortunately, the noise correlation is of known 
functional shape, which is a linear anti-correlation of the 
PAq and RM errors. This allows us to detect this noise 
via a gradient vector-product statistics. 

J d^x p{x) ■ q{x) 



V = V\p,^ 



(7) 



to be regarded as significant co- alignment signals. If we 
use the uncorrected PA maps instead of PAq, we get only 
A = 0.24, 0.35, and 0.36 for the three sources for which 
we have PA maps. These lower values indicates that the 
correlation signal is mostly due to noise in the RM maps, 
which imprints itself to the PAq map during the PA to 
PAq =PA — RM correction. This is verifled by the gra- 
dient vector-product statistics, which leads toV = —0.32, 
—0.36, —0.58, and —0.69 for our set of maps, which clearly 
reveals a preferred anti-correlation of RM and PAq fluc;tii- 
ations. Visual evidence for this can be seen in the upper 
middle panel of Fig. Bl, where diagonal stripes decreas- 
ing from left to right exhibit the presence of such anti- 
correlated noise. As argued above, any real astrophysical 
PAq-RM correlation should be best detectable in A + V, 
for which we get values oi A + V = -0.02, -0.02, 0.11, 
and 0.32. This indicates that at most in one case there 
could be an astrophysical correlation, however, its signifi- 
cance has to be investigated since we cannot always expect 
exact cancellation of V and A for correlated noise. 

In order to suppress the spurious signal, which should 
be mostly located on small spatial scales of the order of the 
observational beam, in the following analysis we smooth 
the RM maps with a Gaussian, leaving hopefully only the 
astrophysical signal carrying large-scale fluctuations.^ 

7. RESULTS 

Applying our statistics, we get for PKS1246-410 A = 
0.07, V = -0.17, A + V = -0.10 (after smoothing to a 4" 
FWHM beam), for Cygnus A we get A ^ 0.07, V = -0.07, 
A + V = 0.00 (FWHM = 2"), for Hydra A ^ = 0.05, 
V = -0.03, A + V ^ 0.02 (FWHM = 1"), and for 3C465 
A = 0.13, V = 0.03, A + V = 0.16 (FWHM = 9.5"). Only 
3C465 shows a marginal signature of source intrinsic co- 
alignment. However, its map size in terms of resolution 
elements is significantly smaller than e.g. the maps of Hy- 
dra A and Cygnus A, therefore a larger statistical variance 
of the alignmcmt measurement is plausible. 

In contrast to the observed maps, our synthetic maps 
are free of correlated noise. This is clearly revealed by 
our statistics, which give for the independent pair of maps 
A = 0.06, V = 0.05, A + V = 0.11, and for the co-aligned 
pair of maps A = 0.89, V = -0.03, A + V = 0.86. This 
illustrates the ability of our statistics to discriminate be- 
tween spurious and astrophysical signals. 

In order to have an estimate of the statistical uncertain- 
ties, we also apply our approach to the swapped RM map 
of PKS 1246-410. We get for the PAoEast-RMwcst com- 
parison A = -0.02, V = -0.19, A + V = -0.21, and for 
the PAq wcst-RMEast comparison A = 0.10, V = —0.30, 
A + V = —0.20. This indicates that the statistical error is 
of the order 6A ~ 0.1 and 5V ~ 0.2 for this dataset, taking 
into account that the swapped maps have less correspond- 
ing pixel pairs compared to the original pair of maps. 

We note that an inspection by eye reveals several sharp 
steps in the RM map of PKS 1246-410, which arc on 
length-scales below the beam size and therefore very likely 

^ These values are anonymized by numerical ordering. As long as the signal is dominated by map imperfections (correlated noise, artifacts) 
we refrain from specifying which map gives which A and V values, since this could be regarded as a ranking of the quality of the maps. We 
think this would be improper, since the observations leading to these maps were driven by other scientific questions than we are investigating. 
Therefore the maps can not be expected to be optimal for our purpose. 

* Smoothing of a circular quantity like a PAq angle is not uniquely defined. 



jd^ximwmi' 

This statistic is insensitive to any existing astrophysical 
PAq-RM correlation, since the latter should produce par- 
allel and anti-parallel pairs of gradient vectors with equal 
frequency, thereby leading to = 0. A map pair without 
any astrophysical signal, which was constructed from a set 
of independent random PA maps, will give V = r ^ —1, 
where — 1 < r < is the correlation coefficient of the noise 
calculated in Appendix B. The statistic F as a test for cor- 
related noise fulfills therefore requirements similar to the 
ones formulated for A, with the only difference that a re- 
quirement like No. 1 is not necessary, since the functional 
shape of noise correlations is known. We expect that corre- 
lated noise leads to a spurious PAq-RM co-alignment signal 
of the order A » \ V\, since for perfectly anti-parallel gradi- 
ents these quantities are identical. Therefore, we propose 
to use the quantity A + V as a, suitable statistic to search 
for a source intrinsic PAq-RM correlation: The spurious 
signal in A caused by the map making process should be 
roughly compensated by the negative value of V, whereas 
any astrophysical co-alignment PAq-RM signal only affects 
A, and not V, since there should be no preference between 
parallel and anti-parallel RM and PAq gradients, leading 
to cancellation in the scalar-product average of V. 

We want to apply our statistics to the RM and PAq 
maps of PKS 1246-410 from Taylor, Fabian & Allen (2002), 
Cygnus A from Dreher et al. (1987) and Perley & Carilli 
(1996), Hydra A from Taylor & Perley (1993), and 3C465 
from Eilek & Owen (2002). Ignoring our considerations 
about correlated noise, an application of our alignment 
statistic reveals co-spatial structures in our set of observed 
maps. We find A — 0.34, 0.56, 0.64, and 0.8,'' which have 
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map making artifacts (see Appendix B). Thus it is uncer- 
tain if this dataset has sufficiently high signal-to-noise to 
make it suitable for a PAq-RM alignment analysis. 

In summary, we do not find any evidence for a significant 
large-scale co-alignmcnt of PAp and RM maps of polarized 
radio sources in galaxy clusters. 

8. SUMMARY 

We investigated if there is evidence for co-aligned struc- 
tures in RM and PAq maps of extended radio sources in 
galaxy clusters as claimed by RB in order to argue for 
source-local RM generating magnetic fields. 

First, we have demonstrated that the null-experiment 
performed by RB was poorly designed for testing the corre- 
lation between PAo and RM in PKS 1246-410. The lack of 
phase coherence in the simulated data resulted in less clus- 
tering in the simulated PAq-RM scatter plots compared to 
the observational scatter plots. Using independent, patchy 
distributions of PAq and RM we show that the correlations 
due to the mutual overlap of the RM and PAq patches pro- 
duce horizontal and vertical chains of clusters as seen in 
the PKS 1246-410 data, whereas co-aligned PAq and RM 
patches produce diagonal stripes. The observed cluster- 
ing therefore favors independent PAq and RM maps as 
expected from foreground intracluster magnetic fields. 

Second, we introduce and apply a novel gradient align- 
ment statistic A. This statistic reveals PAq and RM cor- 
relations regardless of whether they are source intrinsic or 
due to artifacts in the observation or RM map making pro- 
cess. Applying this statistic to a number of radio galaxies 
(PKS 1246-410, Cygnus A, Hydra A, and 3C465) does not 
reveal any significant large-scale co-alignment of PAq and 



RM maps. We find significant small-scale co-alignment 
in all observed map pairs, but we are able to fully iden- 
tify these with map making artifacts by another new suit- 
able statistical test, the gradient vector product statistic 
V . Thus, we introduced two new tools to analyze data of 
Faraday rotation studies of extended radio sources. They 
are powerful in revealing and discriminating observational 
or map making artifacts (by V), and source-intrinsic PAq- 
RM correlations (by A + V), both indicators of potential 
problems for RM based ICM magnetic field estimates. 

Future, sensitive searches for potential, weak source- 
intrinsic PAq-RM correlations with our, or similar statis- 
tics would require observational datasets with a much 
higher signal to noise ratio, and a very well defined obser- 
vational (dirty) beam. We add that such datasets would 
also be crucial for detailed measurements of the magnetic 
power spectra of the ICM, as proposed by Enfilin & Vogt 
(2003). 

We conclude that the observed RM signals of radio 
galaxies embedded in galaxy clusters seems to be domi- 
nated by ~ 1 — 10/LtG ICM magnetic fields in accordance 
with independent evidence. 

We thank Chris Carilli for the data of Cygnus A and 
for his suggestion of exchanging areas of the RM map to 
get uncorrelated maps with conserved higher order statis- 
tics. Data for 3C465 was kindly provided by Jean Eilek 
and Frazer Owen. We acknowledge lively discussion with 
Larry Rudnick. The presentation benefited strongly from 
comments of two anonymous referees, from Larry Rudnick 
and Katherine Blundell, and from Matthias Bartelmann. 



APPENDIX 

PARAMETERS USED FOR THE ARTIFICIAL RM AND PAq MAPS 

The parameters to generate the artificial RM and PAq maps are: TV = 22 for both the RM and the PAq maps. The final 
map size is iVpxi X iVpxi — 50x50pixels. The seed points are distributed by assigning .x™''''* = ([—1.25, 1.25], [—1.25, 1.25]) L, 
where [a, h] denotes a random number, uniformly drawn from the interval from a to b, and L denotes some arbitrary 
length miit so that the final map has coordinates running from —L to +L. Further, RM^ = [—1600, 1600] [0, l]^/^radm~^, 
VRM, = 320[0,l]i/2(cos6'f^,sin6'f^)radm-2L-i, with 6*^^ = [0,27r], and ctrm(x) = [-48, 48] radm-^. Similarly, 
PAi = [-90,90]deg, VPAj = 18 [0, 1]!/^ (cosi9P, sin6'P) degL-\ with 0f^ = [0,2%], and c7pa(x) = [-2.7, 2.7] deg. 
Terms like [0, 1]" lead to non- uniform distributions between and 1, which for < a < 1 are biased towards smaller 
values. 

AVOIDABLE AND UNAVOIDABLE CORRELATED NOISE IN RM AND PAq MAPS 

RM and PAq maps are generated from the same set of radio maps being subject to observational noise. Therefore the 
noise of the RM and PAq maps should be expected to be correlated. In order to understand this correlation one has to 
investigate the map making procedure. 

In order to derive the RM maps, one requires a number n of PA maps at different wavelength Aj, which we denote by 
PAi(x). The first crucial step is to solve the so called 'n-7r' ambiguity, which arises from the fact that the PA is only 
defined modulo tt, whereas the determination of RM requires absolute PA values, since 

PAi(f) = PAo(.f) + RM(f) A? (Bl) 
allows PAj values which deviate more than tt from PAq. Any mistake in solving this ambiguity leads to strong step- 
function like artifacts in RM, and correlated with these, steps in PAq maps. The impact of such steps on any sensitive 
correlation statistics can be disastrous. Fortunately, such n-7r ambiguities can be strongly suppressed by using map-global 
algorithms to assign absolute PAs (Dolag et al., in prep.). 

However, even if the n-w ambiguity is solved, observational noise SPAi in the individual maps, which we assume to be 
independent from map to map, leads to correlated noise in RM and PAq, as we see in the following analysis. Typical RM 
mapping algorithms use a statistic, giving 

RM°^^ = (A2PA - A2 PA)/(A4 - A2^) , and PAg^^ = PA - RM A^ , (B2) 
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where we used Y = X^iLi ^i/*^ denote an average of a quantity Y over the different observed wavelengths. For simplicity 
only, we assume similar noise levels in all maps: = (5PA^) . Here the brackets denote the statistical average and 

should not be confused with the alignment product. One finds for the errors SRM = RM°^'-RM and SPAq = PA^^'-PAq 

(5RM2) = (5PA2)/(n (A4 - A2^)) , (SPAl) = F (^RM^) , and {5RM JPAq) = -A^ (^RM^) , (B3) 
which implies an anti-correlation between PAq and RM noise, with a correlation coefficient of 

r = (^RM (5PAo) /^{SRW^) (SRW^) = VA^ . (B4) 

For the set of frequencies used to derive the RM map of PKS 1246-410 one finds r = —0.89, therefore we expect the noise 
in the RM and PAq maps of PKS 1246-410 to be highly anti-correlated. The scalar product of the gradients, which enter 
the gradient scalar-product statistic V, have expectation values of 

(VRM°^^ • VPAf"'} = {V6RM ■ WPAq) = -{\V6PA\'^)'>^/{n(^ - ^^)). (B5) 

Similarly, one finds 

(|WRM|2) = (|WPA|2)/(n(A4- A2^)) , and {\f6PAof) = {\VSPAf)X^/{n(^->^^)) (B6) 

so that for noise dominated maps there is a strong signal in the gradient scalar-product statistic V = r = — A^/VA^- In 

the case of non-constant astrophysical RM and PAq values, \V\ will be smaller. However, since noise is usually strongest 
on the smallest scales, the gradient of the noise can easily be stronger than the gradient of the astrophysical signal, leading 
to y ~ r. 
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Fig. B1. — Upper left panel: Uncorrected PA-RM scatter plot of PKS 1246-410 (similar to Figure 4 of RB). The underlying maps were 
smoothed by a median weight filter in the same way as RB smoothed the data in order to amplify the clustering. Upper middle panel: 
Similar to the upper left panel, but now corrected FAq-RM scatter plot of PKS 1246-410. Note the appearance of diagonal stripes, which are 
artificial correlations due to errors in the RM maps, which imprint themselves onto the PAq in the Faraday de-rotation step. If they were 
true correlations, no specific direction within the scatter plot should be favored. The decreasing trend from left to right indicates actually 
an anti-correlation, as predicted in Appendix B as a signature of observational noise which propagates through the map maJcing process. 
Upper right panel: PAq-RM scatter plot of PKS 1246-410, but with the RM maps of the Eastern and Western radio lobe exchanged 
as a null-experiment. The clustering remains despite the fact that PAq and RM maps are now independent by construction. Lower left 
panel: Simulated PAq (white lines) and RM (greyscale) maps with independent coherence patches. Lower middle panel: PAq-RM scatter 
plot of the simulated independent PAq and RM maps shown on the left. Lower right panel: PAq-RM scatter plot of simulated maps 
with co-aligned PAq and RM coherence patches. Note the horizontal and vertical patches in the plot of independent PAq and RM. Similar 
structures are clearly visible in the scatter plots of PKS 1246-410 but are absent in the scatter plots of co-aligned PAq and RM. We note 
that these qualitative impressions about the bridges between clusters in the scatter plot are reflected in the results of our two gradient 
statistics introduced in Sects. 5 and 6: The gradient vector product statistic V is sensitive to an imbalance between the diagonally elongated 
features (increasing versus decreasing diagonals which correspond to correlated versus anti-correlated changes in RM and PAq). The gradient 
alignment statistic A measures if both variables RM and PAq change simultaneously as one moves in the maps from patch to patch (and 
therefore in the scatter plot from one cluster to another, as in the lower right panel), or if these changes happen sequentially (e.g. flrst RM 
changes rapidly, and later PAq changes rapidly, as in all upper and in the lower middle panels). 



